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ABSTRACT 

New insights into the formation of interstellar formamide, a species of great relevance 
in prebiotic chemistry, are provided by electronic structure and kinetic calculations 
for the reaction NH 2 + H 2 CO -A NH 2 CHO + H. Contrarily to what previously sug¬ 
gested, this reaction is essentially barrierless and can, therefore, occur under the low 
temperature conditions of intestellar objects thus providing a facile formation route 
of formamide. The rate coefficient parameters for the reaction channel leading to 
NH 2 CHO + H have been calculated to be A = 2.6 x 10 -12 cm 3 s -1 , /? = — 2.1 and 
7 = 26.9 K in the range of temperatures 10-300 K. Including these new kinetic data in 
a refined astrochemical model, we show that the proposed mechanism can well repro¬ 
duce the abundances of formamide observed in two very different interstellar objects: 
the cold envelope of the Sun-like protostar IRAS 16293-2422 and the molecular shock 
L1157-B2. Therefore, the major conclusion of this Letter is that there is no need to 
invoke grain-surface chemistry to explain the presence of formamide provided that its 
precursors, NH 2 and H 2 CO, are available in the gas-phase. 
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1 INTRODUCTION 

Interstellar molecules including a peptide bond -NH- 
C(=0)-, such as formamide (NH 2 CHO) and acetamide, 
(NH 2 COCH 3 ), are particularly interesting for their poten¬ 
tial role in prebiotic chemistry (Saladino et al. 2012 ). For¬ 
mamide was detected for the first time in space towards Sgr 
B2 by Rubin et al. (1971) and later in Orion KL (Turner 
1989). Since then, it has been observed in several massive hot 
cores (Bisschop et al. 2007; Adande et al. 2011), low-mass 
hot corinos (Lopez-Sepulcre et al. (2015), the cold envelope 
and hot corino of IRAS16293-2422 (Kahane et al. 2013), and 
the outflow shock regions L1157-B1 and B2 (Yamaguchi et 
al. 2012; Mendoza et al. 2014). Formamide was also detected 
in the Hale-Bopp comet (Bockelee-Morvan et al. 2000) and, 
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more recently, in comets C /2012 F 6 and C/2013 R1 (Biver 
et al. 2014). Therefore, formamide is present in a large vari¬ 
ety of star-forming environments, as well as in Solar System 
comets, thus supporting the hypothesis that molecules with 
a strong prebiotic potential could have been delivered to 
Earth by comets after being synthesized in prestellar envi¬ 
ronments (e.g. Caselli &; Ceccarelli, 2013). 

In a very recent study by Lopez-Sepulcre et al. 
(2015), formamide has been searched for in ten low- and 
intermediate-mass pre-stellar and protostellar objects as a 
part of the IRAM Large Programme ASAI (Astrochemical 
Surveys At IRAM), which makes use of unbiased broad¬ 
band spectral surveys at millimetre wavelengths. While the 
related species HNCO (isocyanic acid) has been detected in 
all objects, formamide has not been identified in five of them, 
which are the coldest and devoid of hot corinos. According 
to those results, Lopez-Sepulcre et al. (2015) suggested that 
HNCO is formed in the gas phase during the cold stages of 
star formation, while NH 2 CHO is formed from the hydro¬ 
genation of HNCO on the ice mantles of dust grains and re¬ 
mains frozen until the temperature rises enough to cause the 
icy grain mantles to sublimate. Nevertheless, very recent ex- 
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perimental work on the hydrogenation of the frozen HNCO 
seems to dispute the suggestion that HNCO is the parent 
species of NH 2 CHO on ice (Noble et al. 2015). Other het¬ 
erogeneous processes on the icy surface of interstellar grains 
have also been considered (Jones et al. 2011; Walsh et al. 
2012; Garrod et al. 2008). 

The formation routes of NH 2 CHO in the gas phase 
have been only partially investigated. Quan &; Herbst (2007) 
suggested that the radiative association reaction between 
formaldehyde and protonated ammonia followed by dissocia¬ 
tive electron recombination could be a source of formamide. 
However, the model of Quan V Herbst could only produce 
an abundance of formamide of ~ 10“ 15 (with respect to 
H 2 ), significantly lower than the value found by Hollis et 
al. (2006) in Sgr B2. Instead of radiative association, Halfen 
et al. (2011) suggested that the ion-molecule reaction be¬ 
tween formaldehyde and protonated ammonia could lead to 
protonated formamide, which, in turn, forms formamide by 
dissociative recombination according to NH 3 CHO + + e _ —>> 
NH 2 CHO + H. A recent theoretical study of the reaction 
NH 4 + + H 2 CO confuted this suggestion (Redondo et al. 
2014a) because it was verified that all the channels of the 
potential energy surface leading to protonated formamide 
exhibit high-energy barriers. The feasibility of other gas- 
phase ion-molecule reactions that could produce precursors 
of formamide has been recently explored by Redondo et al. 
(2014b) who considered the ion-molecule reactions between 
NH 3 OH+ and NH 2 OH+ with H 2 CO and HCOOH. Also in 
these cases, the presence of high-energy barriers along the 
reaction pathways has been exhibited. 

As for neutral-neutral reactions, the NH 2 + H 2 CO 
—y NH 2 CHO + H seems to be a viable route 
and, indeed, this reaction has been initially consid¬ 
ered in the OSU database (Harada & Herbst 2008; 
http://faculty.virginia.edu/ericherb) with an estimated rate 
coefficient in the gas kinetics range (10 -10 cm 3 s -1 ), with the 
assumption that it is a barrier-less reaction. Later on, how¬ 
ever, Garrod (2013) disregarded such an assumption with 
the following reasons. According to the theoretical study by 
Li V Lu (2002), the more exothermic product channel lead¬ 
ing to NH 3 + HCO is characterized by an entrance barrier 
of 5.89 kcal mol -1 (~ 3000 K) with an estimated rate coef¬ 
ficient of 5.25 x 10 -17 cm 3 s -1 . Li V Lu (2002) investigated 
only the channel leading to NH 3 + HCO. Nonetheless, Gar¬ 
rod (2013), by drawing an analogy with the reaction OH 
+ H 2 CO where the channel leading to HCOOH + H has a 
branching ratio of ~ 2% compared to the dominant H 2 0 + 
HCO channel, concluded that the role of the NH 2 + H 2 CO 
NH 2 CHO + H must be irrelevant. The reaction was there¬ 
fore excluded from the gas-phase network and rather con¬ 
sidered a possible formation route of formamide when in 
ice-assisted chemistry (Garrod 2013; Garrod et al. 2008). 

In this work, we present state of the art quantum me¬ 
chanical characterizations of the characteristic stationary 
points (minima and first order saddle points, also referred to 
as transition states) on the potential energy surface (PES) 
governing the NH 2 + H 2 CO reaction channels (§2 and 3). 
At variance from previous suggestions by Garrod (2013) and 
Garrod et al. (2008), the reaction pathway proceeding with 
the addition of the nitrogen atom of the NH 2 group to the 
carbon atom of formaldehyde has been found to be essen¬ 
tially barrierless. Actually, the entrance channel is charac¬ 


terized by a minimum corresponding to a Van der Waals 
complex followed by a first order saddle point, but both 
stationary points lie below the asymptotic reactants energy 
according to CBSQB3 calculations. We have chosen to treat 
the Van der Waals region adiabatically rather than assume 
that energy randomization takes place at the minimum, as 
such a process is not expected to be efficient. For this rea¬ 
son, we have also employed capture theory as well as Rice- 
Ramsperger-Kassel-Marcus (RRKM) calculations to derive 
the reaction rate coefficient as a function of the temperature 
in the range between 10 and 300 K, that is, the range of rel¬ 
evance for the interstellar objects of interest. The approach 
is the same recently used by us to investigate a bimolec- 
ular formation route of cyanomethanimine, another com¬ 
plex organic molecule with a prebiotic potential (Vazart et 
al. 2015). The derived rate coefficient has also been tested 
in astrochemical models and compared with the formamide 
abundance measured in star formation regions (§4). A final 
section (§5) discusses the implications of the computations 
presented here and the future perspectives. 


2 COMPUTATIONAL DETAILS 
2.1 Electronic calculations 

All calculations have been performed with a development 
version of the Gaussian suite of programs (Frisch et al. 2013). 
Most of the computations were performed with the double¬ 
hybrid B2PLYP functional (Grimme 2006) in conjunction 
with the m-aug-cc-pVTZ basis set (Papajak et al., 2009; 
Dunning 1989) where d functions on hydrogens have been 
removed. Semiempirical dispersion contributions were also 
included by means of the D3BJ model of Grimme (Goerigk 
V Grimme 2011; Grimme et al. 2011). Full geometry op¬ 
timizations have been performed for all molecules checking 
the nature of the obtained structures (minima or first order 
saddle points) by diagonalizing their Hessians. More accu¬ 
rate electronic energies were obtained by the Complete Ba¬ 
sis Set (CBS-QB3) method, which employs a coupled cluster 
ansatz including single and double excitations together with 
a perturbative estimation of triple excitations (CCSD(T)) 
in conjunction with complete basis set extrapolation (Mont¬ 
gomery et al. 2000; Ochterski et al. 1996). 


2.2 Kinetic calculations 


We have performed kinetic calculations using an in house 
code described in previous papers (Leonori et al. 2009, 2013; 
Vazart et al. 2015). The initial bimolecular rate constant 
leading from the reactants to the intermediate is evaluated 
using capture theory calculations, after fitting the energy 
values for the approaching reactants to a 1 /R 6 functional 
form and assuming that each successful capture leads to the 
intermediate. As far as dissociation back to reactants is con¬ 
cerned, we have used a detailed balance argument, whereby 
the unimolecular rate constant for back-dissociation is given 
by the equation 


hback — h C apt 


p{R) 

p(i) 


(i) 


where k cap t is the capture rate constant, p(R) is the den¬ 
sity of states per unit volume for the reactants and p(J) is 



Gas phase formation of formamide 3 


Relative 

Energies 



Figure 1. Proposed reaction path for formamide formation. The 
y axis reports the relative electronic energies (black: CBS-QB3) 
and free energy differences at 0 K (green: B2PLYP-D3/m-aug-cc- 
pVTZ harmonic frequencies). 


the density of states for the intermediate. For the interme¬ 
diate dissociation into formamide + H, where, as opposed 
to the reactants, a well-defined transition state exists, we 
have performed a RRKM calculation. The microcanonical 
rate constant is calculated using the formula 


k(E) 


N{E) 

hp(E) 


( 2 ) 


where N(E) denotes the sum of states in the transition state 
at energy E, p(E) is the reactant density of states at energy 
E and h is Planck’s constant. N(E) is obtained by inte¬ 
grating the relevant density of states up to energy E and 
the rigid rotor/harmonic oscillator model is assumed. Both 
densities of states (reactant and transition state) are appro¬ 
priately symmetrized with respect to the number of identical 
configurations of the reactants and/or transition state. Tun¬ 
neling and quantum reflection have been taken into account 
by computing the tunneling probability for an Eckart bar¬ 
rier having the same negative second derivative at the max¬ 
imum of the pertinent saddle point. After all calculations 
have been performed, the branching ratio between products 
and back-dissociation is determined for each energy and the 
corresponding capture rate constant is multiplied by this ra¬ 
tio to give the rate constant for the formation of formamide. 
Finally, the rate constants are Boltzmann-averaged in order 
to provide the rate constants as a function of temperature. 


3 RESULTS 

3.1 Electronic and vibrational investigations 

The structures of all the molecules were optimized at the 
B2PLYP-D3/m-aug-cc-pVTZ and CBS-QB3 levels of the¬ 
ory. All the precursors, intermediate and products were fully 
characterized as minima on the potential energy surface 
and transition states exhibited a single imaginary frequency 
thanks to vibrational calculations at B2PLYP-D3/m-aug- 
cc-pVTZ level. 

Figure [l] depicts the possible path of the full CH 2 0 + 



Figure 2. Possible radical mechanism for formamide formation. 


NH 2 —> NH 2 CHO formation reaction and the relative elec¬ 
tronic energies of all the involved compounds, obtained at 
the CBS-QB3 level. Zero point energies (ZPE) issuing from 
B2PLYP-D3BJ/m-aug-cc-pVTZ harmonic frequencies were 
added to CBS-QB3 electronic energies to obtain the OK free 
energies (of course identical to enthalpies) also shown in Fig¬ 
ure [l] The proposed pathway includes one intermediate (1) 
and one transition state (2). In the first step of the reac¬ 
tion the NH 2 radical attacks the carbon atom of formalde¬ 
hyde, leading to the H 2 C(NH 2 )0 radical 1 , that is about 
71 kJ/mol more stable than the precursors. This step is fol¬ 
lowed by hydrogen loss from the carbon atom, leading to 
formamide together with H radical, and ruled by the tran¬ 
sition state 2, which exhibits a 67 kJ/mol barrier (reduced 
to 46 kJ/mol if ZPE is included). The products were found 
to be ~ 46 kJ/mol more stable than the reagents. The pro¬ 
posed radicalic mechanism is sketched in Figure [2] 

Alternative paths have been considered, including H mi¬ 
grations in compound 1 (from C and N to O atoms), NH 2 
attack on the oxygen atom of formaldehyde or hydrogen ab¬ 
straction from NH 2 which could have led to NH 3 and HCO. 
However, all these reaction channels are ruled by activation 
energies too high to be overcome in the interstellar medium. 
A plausible tautomeric form of formamide was also inves¬ 
tigated, but once again the corresponding reaction chan¬ 
nel is closed under the typical conditions of the interstellar 
medium. 


3.2 Kinetics 

Figure [3] shows the variation of the computed rate constant 
for formamide formation as a function of temperature up to 
300 K. The calculated rate constants have also been fitted, 
for temperatures from 10 K onwards, to an expression of the 
form: 

k(T) = A x (T/ 300) /3 e“ 7/T (3) 

with A = 2.6x 10 -12 cm 3 s -1 , /3 = — 2.1 and 7 = 26.9 K. The 
rate constant decreases monotonically as the temperature 
increases. 

Even though the activation energy governing H loss is 
marginally higher than the energy of the reactants, this dif¬ 
ference is so small that, even at the lowest energies con¬ 
sidered, tunneling leads to an appreciable rate for hydro¬ 
gen loss. In fact, at low energies (and temperatures), H loss 
largely predominates over back-dissociation. The rates of the 
two processes become equal around 10 K and then, as the en¬ 
ergy increases further, the rate of back-dissociation increases 
much faster than the one for H loss due to the more rapid 
increase of the density of states of the reactants (which in¬ 
clude free rotations) and, as a result, the rate for formamide 
formation drops. 

Figure shows the unimolecular rate constant and 
back-dissociation as a function of energy. 
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Figure 3. Total bimolecular rate constant for formamide forma¬ 
tion as a function of temperature. 



Figure 4. Unimolecular rate constant for dissociation to reac¬ 
tants of initial intermediate as a function of total energy. 

3.3 Astrochemical modelling 

In order to test the impact of the new rates on the pre¬ 
dictions of the formamide abundance, we run a series of 
models with the aim to reproduce the observations in two 
different objects: the cold envelope of the Sun-like proto¬ 
star IRAS 16293-2422 (Jaber et al. 2014) and the molecular 
shock L1157-B2 (Mendoza et al. 2014). We chose these two 
cases because they represent two extreme conditions, with 
the lowest 6 X 10 -12 ) and highest 2 x 10 -8 ) mea¬ 
sured abundance of formamid^] For the model, we used 
the time-dependent Nahoon cod^] (Wakelam et al. 2012). 
It resolves the gas-phase chemical equation as a function of 
the time, using the KID A 2014 chemical network (Wake¬ 
lam et al. 2015), containing 489 species and 7499 reactions. 
We corrected the rate coefficient of the NH 2 + H 2 CO reac- 

1 Note that the abundances here are given with respect to the H 
atoms, while in the original articles they are with respect to H 2 . 

2 http://kida.obs.u-bordeauxl.fr. 



Time (yr) 

Figure 5. Abundances (with respect of H atoms) in the cold en¬ 
velope of IRAS 16293-2422 as a function of time (in years): H 2 CO 
(blue), NH 2 CHO (red) and NH 2 (green). The dashed boxes show 
the relevant measured values. 

tion according the computations presented in the previous 
section. In the following, we will discuss the two cases of 
IRAS 16293-2422 and L1157-B2 separately. 

IRAS16293-2422, cold envelope: IRAS16293-2422 is a well 
studied Sun-like protostar of 22 L® (e.g. Crimier et al. 2010; 
Caux et al. 2011). Specifically, the density of the envelope 
increases going inward with a power law of ~1.5, from a H 2 
density of 2x 10 5 cm -3 at 7000 AU. The temperature also in¬ 
creases going inward from 10 K at the border to 100 K at 85 
AU (Crimier et al. 2010). Kahane et al. (2013) detected the 
formamide using the unbiased spectral survey TIMAS S^] 
(Caux et al. 2011). Successively, Jaber et al. (2014) mod¬ 
eled the formamide line emission to disentangle the contri¬ 
bution from the cold and the warm (hot corino) part of the 
IRAS 16293-2422 envelope and found that the average for¬ 
mamide abundance in the cold envelope is ~ 6 X 10 -12 . For 
our chemical model, we used an average H density of 2 x 10 6 
cm -3 and temperature of 20 K. We run the model assum¬ 
ing that the elemental abundances in Wakelam et al. (2008; 
Table 1, column EA2) are depleted by a factor 10 (C, O 
and N) and 100 (the heavier elements), in agreement with 
previous observations. With these assumptions, the steady 
state abundance of NH 2 CHO, H 2 CO and NH 2 result in ex¬ 
cellent agreement with the observed values, as shown in Fig. 
[5] Note that the implicit assumption is that the cold enve¬ 
lope is formed by material previously present in the placental 
molecular cloud, whose age is > 10 6 yr. 

L1157-B2: L1157-B2 is a molecular shock which is part of 
the molecular outflow system emanating from L1157-mm 
(e.g. Gueth et al. 1996; Lefloch et al. 2012). This system 
has extensively been studied. Briefly, it is composed by two 
major ejection events, at the origin of two molecular shock 
sites: B2, created by the first event, and Bl, created by 
the second one and which is spatially closer to the cen¬ 
tral object. Bl is the best studied of the two molecular 
shocks, as it has (also) been observed in the framework 
of two unbiased spectral surveys programs, CHESS (Cecca- 
relli et al. 2010; http://chess.obs.ujf-grenoble.fr/ ) and ASAI 
(Lefloch, Bachiller et al. in prep.; http://www.oan.es/ASAI), 

3 http://www-laog.obs.ujf-grenoble.fr/heberges/timasss/. 
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which provided the spectral coverage of the 500-2000 GHz 
with Herschel and the 3, 2, and 1mm bands with IRAM- 
30m, respectively. Mendoza et al. (2014) measured the for¬ 
mamide abundance in B2 and Bl, equal to ~ 2 x 10 8 and 
~ 1 x 10 -8 , respectively. Relevant to the simulations re¬ 
ported here, Lefloch et al. (2012) provided constrains on the 
Bl density and temperature, ~ 70 K and ~ 10 4 5 cm -3 re¬ 
spectively. Podio et al. (2014) used the line emission from 
several molecular ions to further constrain the cosmic ray 
ionisation rate in the region, 3 x 10 -16 s -1 , a parameter 
important in chemistry as it influences the timescale of the 
chemical evolution. We adopted the above values also for the 
B2 molecular shock. The model, as in previous works (Po¬ 
dio et al. 2014; Mendoza et al. 2014), consists in a two-step 
procedure: ( 1 ) in the first step, the steady chemical abun¬ 
dances are computed assuming that the elemental abun¬ 
dances in Wakelam et al. (2008; Table 1, column EA2) are 
depleted by a factor 10 (C, O and N) and 100 (the heav¬ 
ier elements); ( 2 ) in the second step, CO, H 2 O, H 2 CO and 
NH^Jare injected in the gas phase, to simulate the passage 
of the shock, and the chemical evolution is followed until a 
few 10 4 yr, even though the age of the Bl shock is evaluated 
to be around 2000 yr and the B 2 a bit older. The CO and 
the H 2 O abundances (with respect to H) are 8 x 10 -5 and 
1 X 10 -4 , respectively. We run several models varying the 
H 2 CO and NH 3 abundances. Figure [ 6 ] shows the abundance 
of NH 2 CHO, H 2 CO and NH 2 obtained assuming H 2 CO and 
NH 3 abundances equal to 3 x 10 6 and 1 x 10 6 , respec¬ 
tively, in agreement with their measurements by Mendoza 
et al. (2014) and Tafalla & Bachiller (1995). The dashed 
areas show the measured values of NH 2 CHO and H 2 CC0 
The predictions obtained using the new values perfectly re¬ 
produce the observations between 1000 and 3000 yr, the 
estimated age of the shock. The predicted NH 2 abundance 
is about 5-9 xlO -8 . 

We also run a few models for the specific case describing 
Bl. We can reproduce the measured formamide abundance 
and predict a NH 2 abundance of ~ 10 -8 . From the Mendoza 
et al. (2014) upper limit, we derive a 3 cr level upper limit on 
the NH 2 abundance of ~ 3 x 10 -9 , which is a factor ~ 3 lower 
than the predicted one. We caution that these upper limits 
suffer of the uncertainty on the excitation temperature, and, 
therefore, further dedicated observations should be carried 
out to better constraint the predictions. 

In summary, the new rate for the reaction H 2 CO and 
NH 2 —» NH 2 CHO allows gas phase reactions fully justify the 
observed formaldehyde abundances, both in the cold and 
warm gas sources, without the need to invoke specific grain- 
surface reactions for that. 


L1551-B2 


c 

o 

~o 

c 

□ 

-Q 

< 



Time (yr) 

Figure 6. Abundances (with respect of H atoms) after the pas¬ 
sage of the shock in L1157-B2 as a function of time (in years): 
H 2 CO (blue), NH 2 CHO (red) and NH 2 (green). The dashed boxes 
show the relevant measured values. 


computations allowed us to propose a reaction path com¬ 
bined to a plausible mechanism concerning this formation. 
Indeed, the first addition step does not involve any barrier 
and can therefore occur in space. RRKM calculations con¬ 
firmed the feasibility of this reaction since once this first 
addition is done, formamide formation largely predominates 
over back-dissociation at low energies. Moreover, we tested 
the impact of these new kinetic data on the prediction of for¬ 
mamide abundance, comparing it with observations in the 
cold envelope of the Sun-like protostar IRAS 16293-2422 and 
the molecular shock L1157-B2 (lowest and highest measured 
abundance of formamide). The results obtained from these 
simulations are in excellent agreement with the observation 
and, therefore, confirmed our computational protocol. 

A more general conclusion is that neutral-neutral gas- 
phase reactions can account for the formation of relatively 
complex organic molecules even under the extreme condi¬ 
tions of ISM. Grain-surface reactions are often called into 
play to explain the formation of complex organic molecules 
because of supposedly missing formation routes in the gas 
phase. However, not all the possible gas-phase routes have 
been actually explored, either in laboratory experiments or 
theoretically, or correctly included in the astrochemical net¬ 
works and models (see, for instance, Balucani et al. 2015). 
Other studies of critical and yet unexplored neutral neu¬ 
tral gas phase reactions performed with the same theoreti¬ 
cal methods employed here can help to fulfill the gap of the 
missing reactions leading to complex organic molecules, es¬ 
pecially when those reactions cannot be easily investigated 
in laboratory experiments. 


4 CONCLUSION AND PERSPECTIVES 

In this paper, we have provided new insights concerning 
the formation of formamide in the interstellar medium. Our 


4 NH 2 is then synthesised by NH 3 , which is the N- element iced 
major recevoir (e.g. Boogert et al. 2015). 

5 The upper limit on the NH 2 in Mendoza et al. (2014) strictly 

applies to Bl, where the NH 2 were observed but not detected. In 
the case of B2, there are no specific observations. 
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